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Abstract 

The quasi-particle model of quark gluon plasma (QGP) is revisited here with a new 
method, different from earlier studies, without the need of temperature dependent bag 
constant as well as other effects such as confinement effects, effective degrees of freedom 
etc. Our model has only two system dependent parameters and surpraisingly good fit 
to lattice results for gluon plasma, 2-flavor and 3-flavor QGP are obtained. The basic 
idea is to evaluate energy density e first from grand partition function of quasi-particle 
QGP and then derive all other thermodynamic functions from s. Quasi-particles are 
assumed to have temperature dependent mass equal to plasma frequency. Energy 
density, pressure and speed of sound at zero chemical potential are evaluated and 
compared with available lattice data. We further extend the model to finite chemical 
potential, without any new parameters, to obtain quark density, quark susceptibility 
etc. and fits very well with the lattice results on 2-flavor QGP. 
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1 Introduction 



The non-ideal behaviour of QGP seen in lattice simulations [1, 2, 3] of QCD, and the elliptical 
flow observed in relativistic heavy ion collisions (RHICs) lead to hot debate on the nature 
of QGP near the critical temperature, T c , for the last few years [4]. Various models such 
as (a) QGP with confinement effects, (b) strongly interacting QGP (sQGP), (c) strongly 
coupled quark gluon plasma (SCQGP), and (d) quasi-particle QGP (qQGP). In model (a), 
confinement effects like bag constant [5], Cornel potential [6] etc. are assumed to give rise 
to non-ideal behaviour of QGP. In (b) [7], effects of bound states of colorless as well as 
colored hadron resonances are assumed to be responsible. In (c) [8, 9], it is assumed that 
QGP near T c is what is called strongly coupled plasma (SCP) and EoS of SCP in QED 
with proper modifications for QCD fits very well lattice data. In (d), QGP is made up of 
quasi-particles with temperature dependent mass [10, 11]. Different varities of qQGP were 
proposed in order to make the theory thermodynamically consistent as well as consistent 
with perturbative and non-perturbative calculations of QCD [12, 13]. 

Here we revisit qQGP with a new method, which involves a less number of parameters, 
to derive EoS of QGP at zero chemical potential (//). Here \x means quark chemical potential 
which is one third of the baryon chemical potential. We start from energy density, rather 
than pressure as done earlier, and derive various thermodynamic properties by fixing the pa- 
rameters of the model. Further we extend the model to include finite chemical potential and 
obtain quark density (n q ), change in pressure (AP) due to finite /i and quark susceptibility 
for 2-flavor QGP. 

2 Phenomenological Model : 

The basic assumption in this model is that quarks and gluons in QGP are not bare quarks 
and gluons, but they are quasi-particles. Because of their collective behaviour in QGP, 
massless partons acquire masses equal to their respective plasma frequencies and become 
quasi-partons. Following the standard procedure of statistical mechanics [14], the grand 
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partition function is given by, 



Q G = £ e -/^-^) 



s,r 



where the sum is over energy states E r and particle number states N s . Now we assume that 
QGP is made up of non-interacting quasi-partons and on taking thermodynamic limit, we 
get, 



(2) 



k=0 



where q is called q-potential and =F f° r bosons and fermions. f3 and z are temperature and 
fugacity respectively, e k is the single particle energy, given by, 



ek = \Jk 2 + m 2 (T,/i) , 

where k is momentum and m 2 is the temperature dependent mass. The main effects of the 
interaction of bare partons, namely collective effects, are taken in (T, ji) dependent mass 
term and treat them as non-interacting quasi-partons. With this assumption, the average 
energy U is given by, 



U =< E r >= 



d i ^ ^ ze k e ^ 



which on taking continum limit and after some algebra, we get, 
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where gf is the degenarcy and equal to g g = 16 for gluons and equal to 2 rif for quarks, n/ is 
the number of flavors. K\ and K 2 are modified Bessel functions of order 1 and 2 respectively. 

3 Thermodynamics (e, P, Cj) of QGP with zero /i: 

Let us first consider the EoS of QGP with zero chemical potential and take z — 1. Hence 
we get the energy density, expressed in terms of e(T) = e/e s , for the quark gluon plasma of 
quasi-partons is 



e(T) = 
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where e s is the Stefan-Boltzman gas limit of QGP, which may be obtained by taking high 
temperature limits of Eq. (4) for gluons and quarks separately and adding them. m g is 
the temperature dependent gluon mass, which is equal to the plasma frequency, i.e, w? g = 
u p = i8~(2^c + n f)- This is obtained from the finite temperature perturbative calculations 
[20, 21]. Note that both gluons and quarks contributes to the gluon mass. However, in 
the case of m q , the temperature dependent mass of quarks, we take m 2 q = ^-rif and no 
contribution from gluons at linear response level, g 2 is related to the two-loop order running 
coupling constant, given by, 

„ (T) = 6n (a - 3(153- 19^) ln(21n(T/A T )) \ 

A ) (33 - 2n/) ln(T/A T ) V (33 - 2n f f ln(T/A T ) )' { ) 

where At is a parameter related to QCD scale parameter. This choice of a s (T) is motivated 

from lattice simulations. With these values of masses with above a s , we can evaluate the 

e(T) from Eq. (4). Note that the only temperature dependence in e(T) comes from a s (T), 

which has the same form as that of lattice simulations [1] with At as a free parameter. 

Pressure can be calculated from the thermodynamic relation e = — P and we get 

T T + Jt T 2 ' [i) 
where P and T are pressure and temperature at some reference points. Note that the stan- 
dard relation of pressure and q-potential is not valid here because of temperature dependent 
We can rederive it using the same procedure from Patria [14], as follows, 
1 
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Therefore, the standard relation of pressure and q-potential is not thermodynamically con- 
sistent in our case. One need extra term arising because of the temperature dependent mass. 



Hence we use Eq. (7) to evaluate pressure. Note that earlier works are all based on modeling 
the extra term so as to fit the lattice data. 

Once we know P and e, c 2 = Qj^- can be evaluated. 

4 Thermodynamics (n q , AP, Xq) of QGP with finite /i: 

Recently, there are lot attempts to simulate QCD with finite 11 on lattice and results are 
reported in Ref. [15, 16] etc. Here we consider the recent results of Allton et. al. [16] for 
2-flavor QGP and try to extend our model to finite \i and explain their results. Again using 
the standard procedure of statistical mechanics, we have, 

<N>=z^-\nQ G , (9) 



which on continum limit and after some algebra, reduces to 
n n 12 ~ . w _, 1 
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Now we modify earlier m 2 q {T) to m 2 q {T, /i) as 



(10) 



m ' q (T,ri) = »—n f (l + -£^) , (11) 

inspired by QCD perturbative calculations [17]. In our case rif — 2 and g 2 is related to 
two-loop order running coupling constant, discussed earlier, but need to be modified to take 
account of finite ii. Following the work of Schneider [18] and Letessier, Rafelski [19], now we 
change T/A T in Eq. (6) as 



hh a &< (12) 

where a is a parameter which is equal to (1.91/2.91) 2 in the calculation of Schneider for 
ii jT < 1 and 1/n 2 in a phenomological model of Letessier and Rafelski. 
From n q , we may obtain other thermodynamic quantities like, 

AP = P(T, ii) - P(T, 0) = f n q dn (13) 

J 

and 

X. = -^U. (14) 
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5 Results 



It is interesting to see that that this simple model very nicely fits lattice data [1, 2, 3] on all 
three systems, namely, gluon plasma, 2-flavor and 3-flavor QGP in the case of zero chemical 
potential. In Fig. 1, we ploted P(T)/T 4 Vs T for pure gauge, 2-flavor and 3-flavor QGP 
along with lattice results. Note that, in the case of flavored QGP, since there is (10% ± 5%) 
uncertainty in P data [3] on taking continum limit with massless quarks [3], we multiply 
the lattice data by the factor 1.1 and is plotted. For each system A T are adjusted so that 
we get a good fit to lattice results. We have fixed Pq from the lattice data at the critical 
temperature T c for each system. Surprisingly good fit is obtained for all systems with At 
equal to 173, 119 and 74 MeV for gluon plasma, 2-flavor and 3-flavor QGP respectively. We 
have taken n/ equal to 0, 2 and 3 respectively for three systems. 

Once P(T) is obtained, then other macroscopic quantities may be derivable from P(T) 
and no other parameters are needed. In Fig. 2, we plotted e/T 4 Vs T/T c for all three systems 
along with lattice results [22] and it fits well without any extra parameters. All the three 
curves looks similar, but shifts to left as flavor content increases. We have taken T c equal to 
275, 175 and 155 MeV respectively for gluon plasma, 2-flavor and 3-flavor QGP. 

In Fig. 3, c 2 is plotted for all three systems, again with lattice results for gluon plasma. 
Reasonably good fit for gluon plasma and our model's out come for the flavored QGP. All 
the three curves have similar behaviour, i.e, sharp rise near T c and then flatten to the value 
close to 1/3. c 2 is larger for larger flavor content. For 2-flavor QGP it is almost coinsides 
with that of gluon plasma. 

In Fig. 4, n q /T 3 is plotted for 2-flavor QGP and compared with recent lattice data with 
out any new parameters. Only the parameter needed is At which is fixed by the results of 
QGP with zero [i earlier. For different values of /x/T c , all curves have similar behaviour and 
very nicely fits with lattice points for T > 1.2 T c . Note that this result is for a = (1.91/2.91) 2 
[18] in Eq. (12). For the other value, a = 1/ir 2 [19], results are not satisfactory, eventhough 
both results coinside for T > 1.5 T c . 

AP/T 4 and Xg/T 2 are plotted in Fig. 5 and 6 and again fits with the lattice results. 
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Note that in the case finite [i case also we multiply the lattice data by a factor 1.1 as we did 
for /i = case. 

Very close to T = T c , i.e., T < 1.2 T c , fits or outcome of our model is not good, especially 
for e, Cg, n q etc. Lattice data also has large error bars very close to T c . However, except for 
small region at T = T c , our results are very good for all regions of T > T c . For T < T c , any 
way our model is not applicable because the system may not be in plasma state. 

6 Conclusions : 

We revisited the quasi-particle picture of QGP with a new method with less parameters. 
In this method we derive the energy density, e(T), first and not P(T) as done earlier, and 
derive P(T), c 2 s etc. QGP is assumed to be consists of non-interacting quasi-partons with 
temperature and chemical potential dependent masses. Interactions between bare partons 
lead to collective effects in QGP and hence they acquire masses equal to plasma frequencies. 
The plasma frequency depends on running coupling constant and we used 2-loop order 
a s (T,n), which is similar to that of lattice simulations, with one parameter related to QCD 
scale parameter. The thermodynamic properties of QGP depends on this parameter as well 
as one integration constant in pressure term which we fix to the lattice point at T c . Both 
the parameters are different for different systems like gluon plasma, 2-flavor and 3-flavor 
QGP etc. Using two system dependent parameter, very good fits to lattice results were 
obtained for energy density, pressure, speed of sound, quark density, quark susceptibility 
etc. No other effects like temperature dependent bag pressure, confinement effects, effective 
degrees of freedom etc. were needed to fit the lattice results. In comparison with other 
models, which generally use more than two system dependent paarameters, here, just with 
two system dependent parameter and using a s (T,/i), inspired by lattice simulation of QCD, 
fits lattice results very nicely. Hence the non-ideal effects seen in lattice simulations of QCD 
may be related to collective behaviour of QGP. Recently, very good fit to lattice results on 
EoS of QGP is also obtained [9] by treating QGP as SCQGP and again the nature of EoS 
is determined by plasma parameter which depends on collective properties of plasma. 
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Figure 1: Plots of P/T 4 as a function of T from our model and lattice results (symbols) for 
pure gauge (lower curve), 2- flavor QGP (middle curve) and 3-flavor QGP (upper curve). 

Figure 2: Plots of e/T A as a function of T/T c from our model and lattice results (symbols) 
for pure gauge (lower curve), 2-flavor QGP (middle curve) and 3-flavor QGP (upper curve). 

Figure 3: Plots of c 2 s as a function of T/T c from our model for pure gauge (lower curve), 
2-flavor QGP (middle curve) and 3-flavor QGP (upper curve) and also with lattice data for 
pure gauge. 

Figure 4: Plots of n q /T 3 as a function of T/T c from our model of 2-flavor QGP and also 
with lattice data (symbols) [16]. 

Figure 5: Plots of AP/T 4 as a function of T/T c from our model of 2-flavor QGP and also 
with lattice data (symbols). 

Figure 6: Plots of Xq/T 2 as a function of T/T c from our model of 2-flavor QGP and also 
with lattice data (symbols). 
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